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ABSTRACT  i  It  le  shown  how  the  method  for  determining  the  transfer  function  of 
a  dynamic  system  with  constant  parameters  developed  by  A.  N.  Sklyarevich 
(8klyarevich,  A.  N.  Operatomyye  metody  v  statistlcheskoy  dlnamlke 
aytomatlchesklkh  slstem.  "NAukan,  M.,  1965)  can  be  extended  to  the  case  when 
the  poles  of  the  transfer  function  of  the  system  and  the  Laplace  transform  of 
the  correlation  function  of  the  Input  signal  do  not  coincide  and  are  of 
order  higher  than  one.  A  procedure  for  determining  the  cross-correlation 
function  of  the  Input  and  output  signals  Is  presented.  This  cross-correlation 
function  Is  used  as  a  basis  for  deriving  a  system  of  algebraic  equations  for 
determining  the  coefficients  of  the  frequency  response  function  and  a 
recurrence  formula  for  their  calculation  Is  presented.  It  is  stressed  that  the 
results  obtained  here  can  be  utilised  in  designing  a  special-purpose  computer 
for  calculating  the  transfer  functions.  Two  examples  illustrate  the 
calculation  procedure.  Orig.  art.  hast  26  formulas.  English  Translation: 

11  pages. 


DETERMINING  THE  TRANSFER  FUNCTION  OF  A  DYNAMIC 
SYSTEM  WITH  CONSTANT  PARAMETERS 

B.  B.  Sodell' 


This  article  examines  a  method  of  finding  the  pulsed  transfer 
function  of  a  linear  dynamic  system  with  constant  parameters;  the 
method  is  based  on  a  comparison  of  the  cross-correlation  function  of 
the  input  and  output  signals  of  the  system  with  the  correlation 
function  of  the  input  signal.  Formulas  are  derived  for  the  case  when 
the  poles  of  the  transform  of  the  input  signal  correlation  function 
and  the  transfer  function  of  the  system  do  not  coincide  and  have  an 
order  higher  than  one. 


In  an  Investigation  of  the  quality  and  accuracy  of  operation 
of  control  systems  during  their  normal  exploitation  the  problem  of 
determining  the  pulsed  transfer  function  g(t)  in  a  time  region  or 
of  the  transfer  function  of  the  system  W(s)  in  the  region  of  a  complex 
variable  can  arise. 

To  find  the  pulsed  transfer  function  g(t)  for  stationary  systems 
with  stationary  random  actions  on  the  input,  different  methods  of 
solving  the  integral  equation 

Kyx  (-)  =-/e  («)  K,  ( •  -u)d"  ( 1 ) 

relative  to  the  function  g(t)  are  used  [1]. 

Let  us  briefly  examine  some  of  them. 
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1.  Calculation  of  the  pulsed  transfer  function  g(t)  by  solving 
the  approximating  system  of  linear  algebraic  equations 

K,A~)  *  i  -  "T^  1  •  ^  2  ) 

Equation  (2)  is  solved  by  any  of  the  known  methods.  Application 
of  the  Oauss-Seidel  iteration  methods  makes  it  possible  to  use 
special-purpose  computers  (synthesizors) . 

2.  Solution  of  equation  (1)  with  a  finite  upper  limit  by  the 

selection  method.  Selecting  a  value  of  g(t)f  we  calculate  the 

Integral  in  the  right  part  of  the  formula  until  the  obtained  value  of 

the  function  coincides  with  the  cross-correlation  function  K  (t) 

yx 

determined  from  experimental  data.  "Control  filters"  are  used  to 
solve  the  problem  by  the  selection  method. 

3.  Determination  of  the  frequency  characteristics  of  the 
system.  The  basis  of  the  frequency  method  is  the  reduction  of 
integral  equation  (1)  to  an  algebraic  equation  with  the  aid  of  a 
Fourier  transform.  As  a  result  we  determine  the  frequency 
characteristics  of  the  system  in  the  form 

r<»— 3 

where  Sx(w),  SyX(J“)  are  the  spectral  and  mutual  spectral  density  of 
the  input  and  output  signals. 

Because  of  the  large  volume  of  computations,  it  is  frequently 
difficult  to  use  these  methods.  Furthermore,  due  to  the  errors  in 
determining  the  correlation  functions  iteration  methods  can  lead  to 
divergent  iterative  processes.  The  application  of  the  frequency 
method  is  connectel  with  the  necessity  of  calculating  the  spectral 
density  functions,  which  is  not  always  convenient. 
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Work  [2]  introduced  a  convenient  method  of  finding  the  transfer 
function  of  a  system  that  does  not  require  calculating  the  cross¬ 
correlation  function  of  the  signals  or  the  corresponding  mutual 
spectral  density.  It  is  based  on  the  concept  of  the  noise  weight 
function  K5(t),  which  is  the  correlation  function  of  the  reaction 
of  the  system  to  an  input  signal  and  the  form  of  a  unit  of  white  noise. 
By  using  the  relationships  that  connect  the  bilateral  Lsplace  transform 
of  the  noise  weight  function  K*(p)  and  the  correlation  functions  of 
the  input  and  the  output  signals  Kjn(p)  and  K*ut(P)  ln  the  form 

KSut(P>  ■  K*(p)Kfn(p)-  W 

and  also  the  connection  of  the  bilateral  Laplace  transform  of  the 
noise  weight  function  and  the  transfer  function  of  the  system  W(p), 
which  has  the  form 

KHp)~W{p)W(-p), 

from  the  known  correlation  functions  of  the  signals  we  can  determine 
the  transfer  function  of  the  system  W(p)  and,  consequently,  the 
pulsed  transfer  function  g(t). 

The  method  is  brought  to  a  specific  result  for  the  case  when 
the  correlation  functions  of  the  signals  of  the  system  have  the  form 

Kln(x)  -  i4#e-B,lT*.(cos p, | x |  +  sin | |),  (5a) 

K  .  (t)  -  V*|t|(cosp,M  +  sin  ft,  h|).  (5b) 

OUv 

In  actuality,  simultaneous  fulfillment  of  equalities  (5a)  and 
(5b)  is  encountered  rather  rarely.  If  the  correlation  functions  of 
the  signals  are  presented  by  different  dependences,  there  can  be 
ambiguity  in  determining  the  transfer  function.  As  A.  N.  Sklyarevich 
has  shown  [3]»  this  ambiguity  can  be  removed  by  including  additional 
information  about  the  system  or  about  its  output  signals,  which 
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significantly  expands  the  possibilities  of  this  method. 

A.  N.  Sklyarevich  proposed  another  method  for  finding  the  transfer 
function  of  the  system  W(s).  It  is  based  on  representing  the  cross- 
correlation  function  of  the  two  nonstationary  processes  Kyx(t,  t^)  of 
the  input  and  output  signals 

K,At,ti)  -  JeV  ~  “)  M  \)du  (6) 

when  t  i  t1,  by  an  expression  of  the  form 

Kyx  (/.  /,)  =  Ig{t-u)  K,  (/,  -  u)  dn,  ( 7 ) 

0 

and  when  t  £  t1,  by  the  expression 

*  t 

(A *1)  —J gif  -  «)  AT,  (/,  —  u)du+fg(t~  u )  Km («  —  /,)  dn.  ( 8 ) 

For  stationary  processes  equations  (7)  and  (8)  can  be  rewritten 
in  another  form  (considering  that  t  -  t1  ■  t): 

K,A t) — j g{u)Kx(x  +  «)<*“  (•  ( 9 ) 

and 

K,Az)=fg(-  +  u)KM(u)du+jg(--u)l<A“)<li'  (OO).  (10) 

When  the  transfer  function  of  the  system  and  the  univariate 
Laplace  transform  of  the  correlation  function  of  the  input  signal 
are  proper  fractions  with  simple  poles  and  Bj  in  the  left  half¬ 
plane  of  the  argument  s  and  when  these  poles  do  not  coincide 
(correlation  resonance  is  not  observed): 

* 

(Ha) 

M 

*(Q  =  Ia/V(/?e«,<°;  RePj< 0),  (llb) 


the  expression  for  K  (t)  when  t  £  0  can  be  written  in  the  form 

y  x 

*„(') -  f  /W*,)'*'1"  + 1 W<V  +  F.i-tyM" • 

f-l  t 


(12) 


where  Fx(s)  is  the  univariate  Laplace  transform  of  the  correlation 
function  Kx(t ) . 

When  the  cross-correlation  function  can  be  presented  in  the  form 

«>0),  (13) 

1-1  / 

after  comparing  the  right  parts  of  formulas  (12)  and  (13)  we  can 
determine  the  expansion  coefficients  of  the  pulsed  transfer  function 

n  (1H) 

rA$i)  +  FA-n 

and,  consequently,  we  can  find  g(t)  in  accordance  with  formula  (lib). 

Let  us  show  how  this  method  can  be  expanded  for  the  case  when 
the  poles  of  the  transforms  of  the  input  signal  correlation  function 
and  the  transfer  function  of  the  system  have  an  order  higher  than 


one. 


Let  us  examine  this  problem:  from  the  known  input  signal 


correlation  function  Kx(t)  and  the  known  pulsed  transfer  function  of 
the  system  g(t),  taking  the  conditions  of  the  physical  possibility 
of  the  system  into  account  [g(t)  =  0  when  t  <  0],  let  us  calculate  the 
cross-correlation  function  of  the  input  and  output  of  the  system 

Kyx(T)- 

The  Laplace  transform  of  the  function  Kx(t)  and  the  transfer 
function  of  the  system  W(s)  can  be  presented  in  the  form 


Z(J=T T 

i— i  xi> 
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In  accordance  with  the  examined  case  we  assume  that  the  poles 


of  the  functions  Fx ( s )  and  W(s)  do  not  coincide: 

■***«!• 

The  correlation  function  Kx(t)  and  the  pulsed  transfer  function  of 
the  system  g(t)  are  connected  with  transforms  (15)  and  (16)  by  the 
inversion  formula 

7  z(s)e*ds, 


where  a  Is  the  abscissa  of  absolute  convergence.  Applying  the 
Inversion  formula  to  expressions  (15)  and  (16),  we  obtain 


L  V  „*-l 


*„(«)- 1  l  Ms-r 

i-i  *--i 


gi « 


(v-l)l' 

—  ■)"  2j  1  B*'-—  (r —  1)! 


*-l  r— I 


(17) 

(18) 


Since 


K  *  ,£*„  <*  —  1  —/)!/! 


(19) 


and 


/ 


(20) 


taking  formulas  (13),  (17)  and  (18)  into  account,  we  obtain  an 

expression  for  the  cross-correlation  function  K  (t): 

y* 


V 1  -f »! _y  Vfl  <V V  V  At  X 


x,i,  .  (r —  i  /)•./! 

(--!)' 


1 


+ 


•  (*,-  ft)-*' 


<*«+;vri/ 
(■:  >0). 


(21) 


Prom  expression  (21)  it  is  evident  that  the  cross-correlation 
function  is  a  combination  of  sums  that  contain  exponential 
corresponding  to  both  the  input  signal  correlation  function  Kx(t) 
(17)  and  the  pulsed  transfer  function  g(t)  (18).  Consequently,  to 
solve  the  reverse  problem,  determining  the  pulsed  transfer  function 
from  experimentally  calculated  correlation  functions,  it  is 
expedient  to  present  the  cross-correlation  function  in  the  form 

*,.<»-  £  l  Z  £  (i>0).  (22 

Calculating  the  coefficients  at  each  power  t  of  the  exponential 
6kT 

components  e  of  the  sought  pulsed  transfer  function,  equating 
them  to  the  coefficients  Mkn  of  the  experimentally  calculated 
cross-correlation  function  (22)  and  introducing  the  designation 


-li  A.(-  D-+'Qr{+/ 

Jail  V— 1 


1 


.MV.  -1, 

j-W’+'J 


(«,+?,)*+/  («, 

we  arrive  at  the  following  system  of  equations: 


**./+!» 


^*A"  (*?  *  )  *  —  Bkp  ®*l  ’• 

Mks-\  (R —  2)  I  -  BkR  an  •+•  BkR~ i  o,i ; 

MkS- 2  (R—  3) !  =  BkR  fl,3  +  BkR-i  On  +  BkR-2 akl ; 


Af„  0!  =  BkRakP  -f-  BkR-i  akR-  x  -f  •••-{-  Bm  o*j  -f*  Bk ,  o,| , 


(23) 


(2*0 


where  R  and  N  are  the  maximal  orders  of  the  pole  6k  in  expressions 
(16)  and  (22).  With  system  (24)  we  can  easily  determine  the  unknown 
coefficients  of  the  pulsed  transfer  function 

B„r=- TP  (25) 

where  a  is  the  determinant  of  system  (24);  Ar  is  the  minor  obtained 
by  replacing  the  r-th  column  with  the  left  part  of  system  (24). 

Prom  system  (24)  we  can  derive  a  recurrent  formula  for 
determining  the  unknown  coefficients 
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(26) 


M^ir — 1)1 —  £ 


In  a  particular  case  we  can  show  that  expression  (1^)  follows  from 
formula  (26)  when  r  ■  v  ■  1  (the  case  of  simple  poles). 

Thus,  we  can  propose  the  following  procedure  for  determining  the 
transfer  function: 

'  using  the  recordings  of  the  signal  realizations  on  the  input 
and  output  of  the  system,  by  any  of  the  known  methods  calculate 
the  correlation  function  Kx(t)  and  the  cross-correlation  function 

Kyx(°  (T  “  °); 

-  approximate  the  graphs  of  the  correlation  functions  by 
analytical  expressions  similar  to  (17)  and  (22); 

-  number  the  powers  of  the  terms  of  the  expansion  of  the  correlation 
function  K  ( t ) ; 

yA 

-  calculate  the  values  of  the  coefficients  from  formula  (23); 

-  determine  the  coefficients  Bkr  from  formula  (25)  or  formula  (26); 

-  present  the  transfer  function  of  the  system  in  the  form  of  (16) 
or,  when  necessary,  transform  it. 

This  method  of  determining  the  transfer  function  is  an  expansion 
of  the  method  proposed  in  work  [3]  for  the  case  when  the  poles  of 
the  transfer  function  of  the  system  W(s)  and  the  transforms  of  the 
input  signal  correlation  function  Fx(s)  have  an  order  higher  than  one. 
This  significantly  Increases  the  possibilities  of  the  given  method 
by  creating  prei equlsites  for  alternating  the  calculation  processes 
in  determining  the  characteristics  of  the  system  during  its  normal 
exploitation. 
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Examples 
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1.  The  input  of  the  object  is  fed  a  random  process  with  the 
correlation  function 

The  behavior  of  the  cross-correlation  function  in  the  region  of 
positive  values  of  the  argument  is  approximated  well  enough  by  the 
dependence 

Kyx (-) -  ~ r*» -jr5"- f-1t (0,25t  +  0. 1 87*>). 

Determine  the  transfer  function  of  the  object. 

The  coefficients  Mkn  (22)  have  the  following  values: 

Af„  *=> — |-i  Mn  —  0,1875;  Mn  =  0,25. 

From  formula  (23)  find  the  values  of  akj+1 

fl"  1  —  2  —  1  -f  2  =~T; 

^“-(-1)°[=T^3+3TT3]“-T; 

,  tJ  1  1  1_  3 

l  l)  [( —  i  — 3)2  .  (—  1  +3)Jj  16* 

Using  formula  (26),  find 

Bu*=  1 ;  flu=*0;  B2\  —  1- 

Thus,  the  transfer  function  has  the  form 

r  (s)  =  ~FTT  +  XsTW ' 

2.  The  input  of  the  system  is  fed  a  random  process  with  the 
correlation  function 

/C,(t)  = 

The  cross-correlation  function  in  the  region  t  £  0  is  approximated 
by  the  dependence 
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Since 


Kl  08  .  .-2124  \ 

85  ,+  3757/ 


+  Hrh  +  ^ 


cos’+(‘lr,+'w 

(4-4-)- 


•j  sin  x  j 


<*< t—e-1 
sin  x  - - j y - ’ 

+  t/ 


cus  ■  - 


we  bring  K  (t)  to  the  form 
yx 


(,) "  {wT'(to + to)  +  eXJ  [(4=7?  +  W+J?  ] 

+ v-t/  (477 + to) + [iTOF + 7TOF  ]} + 

+  |2r-n  +  r^(t|_^.). 


Prom  this  expression  it  is  evident  that 


„  l,l.n  1,1 

"w—  (4-y)4  +  (2  +  jf  '  ”“4—  j+2+j’ 


If  1  1  1  • 

■  (4  +»3  +  (2  —tf  ' 

Afjj  * 

<4+ 

2-/;  *■' 

Determine  the  coefficients 

akj+l 

(23) 

1  ,  1 

1 

1 

0,1  4—  y  +  2+y 

;  a(>  = 

(4— y)* 

(2  +/)” 

1  ,  1 

1 

1 

0,1  4 -fy  '  2— y 

;  fl„  = 

(4  -f  y)2 

(2  — jf  ‘ 

°3'  ~  5  ' 


From  formula  (26)  find  the  unknown  coefficients 

flu  *=  0;  fl|2 *=  1 ;  &]{"—  0;  Bn  —  1 ;  /?3|  =-  2;  ZJ3,  =•  (1 

The  transfer  function  has  the  form 

1,1^  2 
“  (s  +  1  -y?  +  (s  +'  l  +y)‘  +  s  +  2 


2  *l±l__+_2_. 

(i3  +  25  +  2)3  ^  5  4-  2 
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Conclusions 


1.  The  introduced  formulas  allow  us  to  find  the  transfer  function 
of  the  system  with  sufficient  simplicity  for  the  case  when  the  poles 

of  the  transfer  function  and  the  transforms  of  the  input  signal 
correlation  function  have  an  order  higher  than  one. 

2.  The  obtained  results  can  be  used  to  develop  special-purpose 
computers  for  automating  the  process  of  calculating  the  transfer 
function. 

The  author  expresses  his  gratitude  to  Candidate  of  Technical 
Sciences  Ya.  A.  Gel'fandbeyn  for  a  number  of  useful  comments 
expressed  during  discussion  of  the  manuscript. 
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